Invariant Imbedding Equations for Electromagnetic Waves in Stratified Magnetic 
Media: Appfications to One-Dimensional Photonic Crystals 
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We derive the invariant imbedding equations for plane electromagnetic waves propagating in 
stratified magnetic media, where both dielectric and magnetic permeabilities vary in one spatial 
direction in an arbitrary manner. These equations allow us to obtain the reflection and transmission 
coeflicients of the waves and the field amplitudes inside the media exactly for any polarization and 
incident angle of the incoming wave by solving an initial value problem of a small number of ordinary 
differential equations. We apply our results to one-dimensional photonic crystals, where the periodic 
variations of both dielectric and magnetic permeabilities create photonic band gaps in the frequency 
spectrum. 



The propagation of electromagnetic waves in inhomo- 
geneous media is an important topic in various branches 
of physics, including optics, plasma physics, astrophysics, 
and condensed matter physics |2| ||, Great atten- 
tion has been paid to the cases where the inhomogene- 
ity is random or periodic @, |, |, |, |7[ |, |, 0. The 
problem of electromagnetic wave propagation in random 
dielectric media is analogous to the Anderson localiza- 
tion problem of noninteracting electrons in a random po- 
tential [|, |ll|, In recent years, much effort has 
been made to demonstrate the localization of light ex- 
perimentally |l4[ |l5| and to use the effect for developing 
the so-called random laser |l^. Another important 
area of research has been created from an analogy with 
the electronic band structure problem in condensed mat- 
ter physics. In media where the dielectric permeability 
varies periodically in space, a series of photonic band 
gaps can appear in their frequency spectra jlj, 
Because of the promising possibility that a number of de- 
vices based on this phenomenon can be used in various 
optoelectronic applications, research activity in photonic 
band gap structures has increased explosively in the last 
decade ll|,|20). 

A number of analytical and numerical techniques have 
been developed for analyzing the wave propagation in 
inhomogeneous media |, |, | | ||, §, p p. In 
the most general cases where the inhomogeneity is three- 
dimensional, the exact solution of the problem is ex- 
tremely hard to obtain, and one is forced to use some 
approximation scheme. In the cases where the inho- 
mogeneity is one-dimensional, however, there exist sev- 
eral theoretical methods that allow (at least numerically) 
exact solutions of the problem in some simple situa- 
tions. The invariant imbedding method is a powerful 
and most versatile tool for handling one-dimensional sit- 
uations 0, |2^, |2^. It can be used to obtain 
exact solutions for the reflection and transmission coef- 
ficients of incoming waves and the electric and magnetic 



field amplitudes inside arbitrarily inhomogeneous media. 
When the inhomogeneity is random, it can be used to 
obtain the exact disorder-averaged reflection and trans- 
mission coefficients and field amplitudes [ p7[ . The in- 
variant imbedding method can also be used to study the 
wave propagation in nonlinear dielectric medi a ||2^ , 
and the propagation of several coupled waves ]24| in an 
exact manner. 

In the present letter, we apply the invariant imbedding 
method to cases where both the dielectric and magnetic 
permeabilities vary in one direction in space in an arbi- 
trary manner and derive a set of first-order ordinary dif- 
ferential equations called the invariant imbedding equa- 
tions. To the best of our knowledge, these equations 
have never been derived before. We solve the initial 
value problem of the invariant imbedding equations nu- 
merically to obtain the exact reflection and transmission 
coefficients and field amplitudes. We apply our results 
to the cases where one-dimensional periodic variations 
of both the dielectric and magnetic permeabilities create 
photonic band gaps in the frequency spectrum. This sit- 
uation can be realized in one-dimensional magnetic pho- 
tonic crystals made of alternating layers of dielectric and 
ferrite films 123, 



3q, 31 



We are interested in the propagation of a plane, 
monochromatic, and linearly-polarized electromagnetic 
wave of angular frequency to and vacuum wave number 
ko — Lj/c, where c is the speed of light in vacuum. The 
wave is assumed to be incident from a homogeneous re- 
gion on a layered or stratified medium, where the di- 
electric permeability e and the magnetic permeability fi 
(and therefore the refractive index n = y^eT* cind the 
wave impedance Z — y/ fj,/e) vary only in one direction 
in space. We take this direction as the z-axis and as- 
sume the inhomogeneous, but isotropic, medium lies in 
< z < L. Without loss of generality, we assume that 
the wave propagates in the a;2;-plane. Since the medium 
is uniform in the x-direction, the dependence on x can be 
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taken as being through a factor e"^^, with q being a con- 
stant. When g = 0, the wave passes through the medium 
normally. If g ^ 0, the wave is said to pass obliquely. 

For g ^ 0, we have to distinguish two independent 
cases of polarization. In the first case, the electric field 
vector is perpendicular to the xz-plane. Then one can 
show easily from Maxwell's equations that the complex 
amplitude of the electric field, E — E{z), satisfies 



(PE 1 dfi dE 
dz^ fi{z) dz dz 



+ [kle{z)ti(z) ~q^]E = 0. (1) 



This type of wave is known as an s (or TE) wave. In 
the other case, the magnetic field vector is perpendicu- 
lar to the a;z-plane. Then the complex amplitude of the 
magnetic field, H = H{z), satisfies 



d^H 



1 de dH 



dz"^ e(z) dz dz 



+ [kle{z)^,{z) -q^]H = 0. (2) 



This type of wave is called p (or TM) wave. We note 
that there is a simple symmetry between Eqs. (^ and 
(^ : Eq. (H) is obtained from Eq. (|]) by exchanging e and 
/i and replacing E by H. From now on, we will consider 
only the s wave. The results for the p wave follow from 
those for the s wave in a trivial manner. 

We assume that the wave is incident from the region 
where z > L and is transmitted to the region where 
z < 0. The dielectric and magnetic permeabilities are 
assumed to be given by 




iei{z) 



if z > L 

if < z < L , 

if z < 

if z > L 
if < z < i 
if z < 



(3) 



where ei, Mi; ^i^d ^2 are real constants and e/j:(z), 
e/(z), fj,ji(z), and fJ.i{z) are arbitrary real functions of z. 
When 9 is defined as the angle of incidence, the constant 
q in Eqs. (^ and (^) is equal to y'ci/iifco sin 6*. 

We consider a plane wave of unit magnitude E{x, z) = 
_E(z)e*^^ = gV{L-z)+iqx ^ where p = ^ei/iifco cos 0, inci- 
dent on the medium from the right. The quantities of 
main interest are the complex reflection and transmis- 
sion coefficients, r — r{L) and t ~ t{L), defined by the 
wave functions outside the medium: 



E{x, z) = 



gip(L-2)+iija; _j_ j,^2^^gip(2-L)+ig£c^ Z>L 

i(i)e-*?''^+*?^, Z < 



(4) 

where p' = ^£2/^2^0 cos 6' and 9' is the angle that outgo- 
ing waves make with the negative z-axis. 

We generalize the invariant imbedding method devel- 
oped in Ref. 23 to the cases where both e and fi are 
inhomogeneous. Let us consider the E field in Eq. (|l|) as 



a function of both z and L: E = E{z;L). Starting from 
the crucial observation that the boundary value problem 
of the wave equation, Eq. (|^), can be transformed into 
an integral equation 



E{z;L)^g{z;L) + 



ip 



dz' g (z; z') 



e(z') m(^') , 



ei 



/ii 



e{z') 



g{z;z') = exp 
we derive 



ip sgn(z — z') 



dz 



^i{z') 

Ml 



dE{z;L) 
dL 



a{L)E{z;L), 



E{z';L), 
, (5) 

(6) 



where 



a{L) 



ip- 



^i{L) ip ( e{L) ^{L) 



Ml 2 1^ ei 
2 r,(£) 



Ml 



ei 



M(i) 



Ml 

E{L;L) 



(7) 



Using this equation, we obtain exact differential equa- 
tions satisfied by r and t: 

dr{L) nR{L) +ifii{L) 

— 2iJeiiiikQ cosW riL) 

Ml 

eniL) + iei{L)^ _ f^{L) + ifiiiL) 
Ml 



dL 

i 



jVeiMi^o cos 6* 

_^ eR{L) + iei{L) 
ei 

x[l + r(L)]\ 
dt{L) 



ei 



tan^^ 



dL 



i^ei/iifco cos 9 



fJ,B.{L) + ini{L) 
Ijlr{L) + ifii{L) 



Ml 

tan^e 



Ml 



t{L) 



2 V^^iMi^o cos 6* 

^ eR,{L)+iej{L) 
ei 

X [f + r(L)]i(L). 



eR{L) + iei{L) uniL) + ii.ii{L) 



ei 



tan^ 9 - 



Ml 



/ifl;(L) -I- iiJ,i{L) 



Ml 



tan 



(8) 



These invariant imbedding equations are supplemented 
with the initial conditions for r and t, which are obtained 
using the well-known Fresnel formulas: 



r(0) 

m 



M2>/eiMi cos^ — fii V 62^2 — sin^ 9 
M2V^iMi cos6' + Ml V e2M2 - eiMi sin 9 

2m2\/^iMi cos6' 
M2 V^iMi cos -I- Ml V e2M2 - eiMi sin^ 9 



(9) 



For given values of ei, 62, Mii M2j (or lo), and 6 and 
for arbitrary functions ei^(L), ei{L), fj,R{L), and fii{L), 
we solve the nonlinear ordinary differential equations in 
Eq. (^) numerically, using the initial conditions in Eq. (|^), 
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FIG. 1: Transmittance spectrum of a 100-period array of al- 
ternating layers with ea = 4, fia = 2, et — 2, Hb — 1, and 
da = di = A/2. 



FIG. 2: Transmittance spectrum of a 100-period array of al- 
ternating layers with ea — 4, fia ~ 1, — 2, fit ~ 2, and 
da = db = A/2. 



and obtain the reflection and transmission coefficients r 
and t as functions of L. The reflectivity TZ and the trans- 
missivity T are given by 



T = 



\t\' 



(10) 



If the permeabilities have no imaginary parts, the quan- 
tities TZ and T satisfy the law of conservation of energy, 
n + T =1. 

The invariant imbedding method can also be used in 
calculating the fleld amplitude E{z) inside the inhomo- 
geneous medium. We rewrite Eq. as 



dE{z-l) 



dl 



iy^eifiiko cosl 



11r{1) + ifii{l) 



Ail 



€R{l)+i€l{l) IJ:r{1) +ifXl{l) 



ei 

eR{l) + iei{l) 
ei 



Ml 



tan^ 



Ml 



tan^ 61 



X [l + r{l)]\E{z;l). 



(11) 



For a given z {0 < z < L), the fleld amplitude E{z; L) is 
obtained by integrating this equation from I — z to I = L 
using the initial condition E{z; z) = 1 + r{z). 

We demonstrate the validity and utility of our invari- 
ant imbedding equations by applying them to simple one- 
dimensional photonic crystals, where both e and fi vary 
periodically in one direction in space. More speciflcally, 
we consider periodic arrays of alternating layers, char- 
acterized by the dielectric permeabilities ea and e;,, the 
magnetic permeabilities fj,a and fib^ and the thicknesses 
da and db- We assume that the imaginary parts of the 



permeabilities are zero and the ambient medium is vac- 
uum. In the limit where the length of the system, L, 
diverges to infinity, one can find the exact positions of 
the photonic band gaps by computing the roots of the 
analytical expression js^ 

cos /3a cos /3h - i ( ^ -I- ^] sin/3aSin/3{, -(-1 = 0, (12) 
2. \nb ha 



where — 



a, h) and 



^ [ ni^ — aiv? 9 / (riiZi) for the s wave 
I Zi\/ ni^ — sin^ /rii for the p wave 

P, = kod,^Jn,^~sm^e. (13) 

In Fig. 1, we plot the transmittance T (= |tp) of s and 
p waves impinging on a 100-period array with ea — 4, 
/ia = 2, Efc = 2, = 1, and da = db as a function of 
fcpA where A — da + db- In this case, the refractive index 
alternates periodically between Ua = 2-\/2 and Ub — V2, 
whereas the wave impedance is uniform. We find that 
there is no photonic band gap when the wave is normally 
incident. For obliquely incident waves, small gaps are 
created. Furthermore, we observe that the band gaps for 
the s wave agree precisely with those for the p wave in 
their sizes and positions. When B — 45° and L oo, 
we find from Eq. (^2|) that the lowest band gap exists 
for 1.54 < fcoA < 1.63, which is fully consistent with our 
numerical result obtained using the invariant imbedding 
method. 

In Fig. 2, we plot T for the case of a 100-period array 
with ea = 4:, fia = 1, eb ^ 2, fib = 2, and da = db. 
In this case, the wave impedance alternates periodically 
between Za = 1/2 and Zf, = 1, while the refractive index 
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is uniform. We find that large band gaps are formed for 
all incident angles. Similarly to the previous case, the 
band gaps for the s wave agree with those for the p wave 
in their sizes and positions. When 6 = 0° and 45°, the 
first gaps are predicted to exist for 1.23 < k^A < 1.91 
and for 1.32 < koA < 2.04, respectively, in agreement 
with our result. 

From the two examples discussed above, we conclude 
that the wave impedance, not the refractive index, is a 
crucial parameter mainly responsible for the creation of 
photonic band gaps, as has been extensively discussed 
in Refs. 29 to 31. The equality of the s and p wave 
gaps is a special feature that works only when either the 
refractive index or the wave impedance is uniform. It is 
easy to verify this by examining Eq. (|l2|). 

It is straightforward to apply our equations to more 
realistic cases where the imaginary parts of the perme- 
abilities are nonzero. They can also be applied to cases 
where the parameters contain randomly-varying compo- 
nents. This versatility of the invariant imbedding method 
makes it a very useful and convenient tool in designing 
practical optical devices for various application purposes. 
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